Interpolation, projection and hierarchical bases in 
discontinuous Galerkin methods 



Lutz Angermann* and Christian Henke^ 
January 12, 2013 



Abstract: 
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directed to applications to discrete conservation laws. 
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1 Introduction 

The topics covered in this paper belong to the fundamentals of the analysis of discontinuous Galerkin 
methods for partial differential equations. On the one hand, we have collected and reproved results 
from areas that are important for any study of finite element methods, such as the piecewise 
polynomial approximation in Sobolev spaces, quadrature formulas and inverse inequalities (Sections 
[2| [3| |4) [§. On the other hand, we have directed our attention to facts that are specifically related 
to particular techniques such as certain relations between lumping and quadrature effects or the 
investigation of fluctuation operators and shock-capturing terms by means of a hierarchical basis 
approach (Sections |4| [5]). A major concern of our study was to trace the dependence of the constants 
on both the mesh width and the local polynomial degree. 

The paper is organized as follows. After a brief introduction, which introduces the basic notation, 
we investigate polynomial approximations by means of tensor product elements on affine partitions 
in Section [3} This includes estimates of the reference transformations, which we prove by the help 
of a general chain rule. In this way we get a certain insight into the structure of the occurring 
constants. After that we prove error estimates for the Lagrange interpolation and the -L2-projection 
both with respect to the elements and with respect to the element edges in the scale of Sobolev- 
Slobodecki spaces. 

In Section [4j we prepare some important notions such as quadrature formulas, lumping operators 
and discrete L2-projections for later purposes. In particular, we point out the importance of an 
suitable choice of the quadrature points for optimal (w.r.t. the local polynomial degree) error 
estimates of the Lagrange interpolation. 
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In the following section we investigate the projection and interpolation errors for Gauss-Lobatto 
nodes. For this purpose we extend the concept of the hierarchical modal basis to the so-called 
embedded hierarchical nodal basis and we prove error estimates for the Lagrange interpolation and 
for the discrete -L2-projection which are optimal on the elements and almost optimal on the element 
edges. 

As a natural complement to the (direct) estimates from the previous sections, we present in Section 
[6] inverse inequalities that are based on generalizations of the Nikolski and Markov inequalities. 

2 Basic notation and definitions 

Let Q C M'^, d E N, be a bounded polyhedral domain with a Lipschitzian boundary (see, e.g., 
|AF03l Def. 4.9]). O is subdivided by partitions T (in the sense of [EG041 Def. 1.49]) consisting 
of tensor product elements T (closed as subsets of M*^) with diameter Ht '■= maXx.y^T \\x — yWfi- 
Here and in what follows the symbol || • \\iP, p £ [l,oo], denotes the usual £P-norm of (finite) real 
sequences. Furthermore the maximal width of a partition T is defined hy h := maxj/iy : T S T}- 
To indicate that a particular partition has the maximal width h, we will write Th- 
in this paper, the standard definition of finite elements {T, Pt,T,t} is used, see e.g. |EG04| Def. 
1.23]. The finite element space is defined by 

Wh := {w G L\n) : w\t e Pt VT G %} , 

where 

PtCW^'°^{T) yreTh for some / > 0. 
Because of the last requirement, Wh is a subspace of 

W'^PiTh) := {w G L\n) : w\t G W'^^iT) VT G T^} , p G [l,oo]. 

A particular finite element {T, is generated by means of a reference element {T, P, S}, 

where the geometric reference element T is mapped onto the geometric element T by a C^- 
diffeomorphism Ft : T ^ T. 

In the case of a Lagrange finite element {T, P, S} (in the sense of |EG04| Def. 1.27]) with the node 
set J\f := {xi, . . . , x^k }, r := and the linear forms 

dof 

ai{v) := v{xi), l<i< n^^j, \/v G P, 

the definitions 

Pj. ■= {voF^^ : u G P} (1) 

and 

ai{v) = ai{v o Ft) = ai{v), v{x) = (v o Fj:^){x) = v{x), x = Ft{x), 1 < « < n^of, 

are used. The nodal basis of Pt is obtained by an analogous transformation of the reference shape 
functions 

{'^i)--->'^n* } with (pi{xj) = 5ij, l<i,j<n\oi- 

dot 

DEFINITION 1 A partition Th is called affine if the mapping Ft is ajfine for all T G Th, i-c. if 

Ft{x) = Jtx + 6r with Jt G M'^-'^, det(JT) / 0, 6t e 

In addition, the following properties of Th are important. 
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DEFINITION 2 (locally quasiuniform, shape regular) A family of affine partitions {Th} h>o 
is called locally quasiuniform if there exists a constant ctq > such that 

V/i > VT G Tft : aT:= — < ctq, (2) 

QT 

where qt denotes the diameter of the largest hall contained in T. 

DEFINITION 3 (quasiuniform) A family of partitions {Th}h>o is called quasiuniform if it is 
locally quasiuniform and if a constant Cqu > exists with 

V/i > VT G 7^ : hT> Cquh. (3) 



3 Polynomial approximation using tensor product elements 



In what follows we will investigate thoroughly the approximation of functions on affine partitions. 
We set f := I := [-1, 1] and P := Qk{f) with 

Q^(f ) := span {x"}, £ G f , A; G No := N U {0}. 
A function ipa £ {'^ij • • • > <^„fe " ^ can be written as a product of univariate Lagrange 

dof 

polynomials. Namely, denote by {ipQ, ipi, . . . , (^^} a basis of the space of univariate polynomials of 
maximum degree k. Then, for any multiindex q G Ng with ||a||^oo < k, 

The space Qk{T) is defined according to ([T]). Similarly, based on the definition 

Pfc(f) := span x G f. A; G Nq, 

"eN;^,||a||^i<fc 

the space ^k{T) of polynomials of maximum degree k can be introduced. 

As a consequence of the affine structure of the transformation between the reference element T and 
the element T we can prove the following estimates. 

LEMMA 1 For I > and 1 < p < oo, l/oo := 0, there exists a constant Ci^d ^ 1 such that, for 
T ^ Th, Th o^ffine, and w G W^''^(T), w = w o Ft, the following estimates hold: 



\w 



l,p,T 



< CiJJTlUdetiJr) 



-111/ 



\det{JT)\^/P\w 



'W\l,p,T, 



(4) 
(5) 



\w\l,p,T < Ci^dWJx 
where Ci^d depends only on I and d. In particular, Cg^rf = 1. 

Proof: By Faa di Bruno's formula (see, e.g. |Joh02] ). for w = w o Ft, \a\ =1, P:=<\{0} and 



A:= 

we have the representation 

d°'(woFT){x) -- 



: P ^ a(7) = /3 and |a(7)|7 = « 

7eP 7eP 



«! E|/3|<|„| {d^n,){FT{x)) EaeA nU ^kv- 



(d''FT){x) 



7! 

{d,FT) 



a(7) 



a{ej) 
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because the conditions I = \a\ = IJ^-yeN"^ \^(.l)h\ = \Yl'j=iWi^j)\^j\ = J2'j=iWi^j)\ ^^"^ l/^l = 
I Yl'j=i — Yl'j=i imply that |a| = = / if "empty" sums are neglected. The absolute 

value of the left-hand side can be estimated as 



mwoFT){x)\ < alZ\^^^,\{d^w){FT{x))\EaeAUU4j 
< alE\,\=i \{d^w){FT{x))\ Ea^AllU ^ 



= c„,,E|/3hI(5M(^t(^))|||Jt||Lx 
< c^4Em=i\id^yj)iFT{m\\JT\\[2. 



(e,)! 



Jt\L 



For p < oo, using Holder's inequahty for sums and observing that X]|^|=z 1 — {'^~^\~^)' 
reference element we have that 



0,p,T — 



d + l-l 

I 



p 

o,p,r' 



Applying the substitution rule, we get 



I 



C^\JT\\t\det{JT)nw\l^^T- 



Summing up w.r.t. a, the estimate 
\w\P ^ = Z 



'l,p,T 



^ ^\a\=l 



0,p,T 

/d+«-l\(P-l)^P II T ||/P| 1 , / T X|-l| |P 



< ( max {Ca4) 

Ja\=l 



CmJT\\%\det{JT)nw\l^T 



)'^||JT||^^|det(JT)|-V|f,,,^ 



proves the statement with Ci^d { ) t^Slx {Ca,d) ■ For p = oo we see that 

\a\=l 



W 



Loo,T 



max \\o w\ 



\a\=l 



0,oo,T 



< max(C7,,,) ||Jr||^.|| E|/3|=/ K^M llooot 

\a\=l ' ' ' ' 

= ma^{Ca,d) WJrim E|/3|=/ (^)\ IIo,oo,t 

\a\=l ' ' 

< max(C7,.^)||JT||^2C'+l"')||max|a^u;| ||o,oo,T 

\a\=l ■ \I3\=1 

= Q,d||'/T||i2 max||a''t(;||o,oo,T 
lit , ^ 



Since Ft : T — )• T is bijective, the second estimate follows obviously. -4 
LEMMA 2 The following estimates are valid: 

|det(JT)| = ^, \\JTy<— and ||J^i||,. < ^. 

|T| pf PT 



(6) 



Proof: The first relation is classical. The proof of the inequalities is easy, see, e.g., |Cia91^ Thm. 
15.2]. < 
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COROLLARY 1 Given a locally quasiuniform family {Th}h>o of affine partitions with the refer- 
ence element T = I'^. Then: 

II Ml < II r-i|l <: ^^^0 ^7^ 
\\JT\\e2 < — , \\Jj. \\i2 < — [7) 

Z tlx 

and 

d 

det(JT) = n^'^' ^ = II^tII,". < 
|r| = |f ||det(JT)| < /iT, 
where Aj, 1 < i < d, are t/ie eigenvalues of the matrix J^Jt- 

Now we are ready to define the Lagrange interpolation operator as follows: 

i=l 

We mention that the results from this section on the Lagrange interpolation operator do not 
impose any conditions w.r.t. the location of the nodes. Later we will formulate statements about 
the interpolation operator which use the special properties of Gauss-Lobatto quadrature points. 
Simple calculations show that 

I^v = V yveQkiT), (9) 
||/Ml/,p,T < C'||i;||o,oo,T, />0, [l,oo], (10) 

where C = C{T,{ifi}^^'f,n^^f,l,p) > 0. Furthermore we have the following error estimates. 



LEMMA 3 Let T be an element of an affine partition Th such that the corresponding family of 

partitions {Th}h>o locally quasiuniform. Assume that l</</c + l, / gN, [l,oo] and Ip > d. 

jk 

such that 



Then, for the Lagrange interpolation operator there exist constants C > independent of Ht 



\v-4v\r,p,T < C/i^-''blz,p,T, 0<r<l, (11) 
\v - ljlv\r,p,E < Ch^-^/P-'lvli^p^T, 0<r<l-l/p, (12) 

for all V E W^'^{T), where E denotes a face ofT. 

The proof relies on the following interpolation inequality. 

LEMMA 4 (interpolation inequality) Let lo,li E Nq, Iq ^ h, 1 < p < oo and G be a bounded 
domain with a Lipschitzian boundary. Define, for < 9 < 1, 1$ := {1 — 9)Iq + 6li. Then 

Vu E VF''"^(G) n W'^'P{G) : \\u\\i,,p,G < CMUJuWI^^g ■ (13) 

Proof of Lemma |4| We first mention that the statement for the case G = M'^ is a consequence of 
|Tri78[ 1.3.3 (g)] and |BL76[ Def. 6.2.2, Thm. 6.2.3, Thm. 6.2.4 and Thm. 6.4.5 (3), (4)]. If G is 
a bounded domain with a Lipschitzian boundary, then there exists a total extension operator Eq 
(see |AFn3[ Thm. 5.24] or [S^eTOl Ch. 6, Thm. 5]) such that 

\\u\\io,P,G = \\Egu\\i,,p,g (by [SEOSl (i) in (5.17)]) 

< \\EGu\\ig,p,Rd 

< C\\Egu\\]-;^^,\\Egu\\1^^^, 

< CWuWj-lJuWlj,^^ (by ISFQSI (ii) in (5.17)]). ^ 
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REMARK 1 Using further results on Stein's extension operator \Kal8^ p. 185 and Thm. 1], the 
interpolation inequality (13) can he extended to the parameter set lo,h £ 7^ ^ij 1 ^ POiPi ^ 

oo, as follows, where : 



1 — U I O 
PO Pi' 



Vn G Ty'o'^^'HG) n Ty'i'Pi(G) : hlb„p«,G < C\\u\\l 



lo,po,G\\^\\h,Pi,G ■ 



Proof of Lemma [sj By the triangle inequality, (10), and the embedding theorem |AF03l Thm. 
4.12, Part II], for 1 < I < k + I, r £ Nq, r < I, Ip > d we have that 



\v — itv 



< \\v\ 



< \\v\ 



+ C\\v\ 



< C\\v\\i,p^f G VF''P(f). 
Furthermore, by Q and Pfc(f ) C Qfc(r), 



0,oo,T 



(14) 



inf_ \\{Id-l'f^){v+: 

P&k{T) 



W,p,T 



< inf . \\{Id-I^)iv+y,, , 



where Id denotes the identity operator. From ( 14 ) we see that 



\v-I^v 



<C inf + ^<C7|{)|^ ^, 



h^\r,p,T — ~ II - ■ r- ii( p 



(15) 



where the last estimate is a consequence of the Deny-Lions lemma ( |Cia9H Thm. 14.1]). Note that 
the constant C depends on the parameters of the reference element. The application of Lemma [T] 
results in an estimate on the element T: 

\V - lj:v\r,p,T < C\\J^%,\det{JT)Mv-lj;vl^^^f 



< C\\J^^\\'^,\detiJT)\^/P\v\ 



Up,T 



< C(||JT||,2||J^i||,2)n|JTH'|nH;,p,r 

< C 



where we have used the condition ([2]) in the last step. 
In the case r G ]R_|_ \ No, we apply the interpolation inequality (13) with Iq := [r 



li := \r~\ := min {z} and 9 := r — [rj: 

z>r 



max z, 

z£Z,{z}<r 



< Ch'f'-lvli^p^T, \r] < I. 



The proof of the second estimate (12) runs similarly. 

For 1 < / < A; + 1, p G (1, oo), and any face E of the reference element T with E = FtE we have, 
by the special trace theorem for the faces of T (see |Nec67t Thm. 2.5.4]) and (14), that 



\v-l'^v\ 



< cwv-nm,^^ < c\\v\ 



Vt) G W^^P{f). 



'h^\\l-l/p,p,E - ^W'^ ^h^Wl^p^T - ^ll''ll«,p,r 

Furthermore, |Nec67t Lemma 2.5.4] implies that, for any r G [0, Z — 1/p), 

\\v - /Ml,,,,^ <C\\v- lti\\i-i/p,p,E V* G W^'^if) 



(16) 



Combining this estimate with (16) we arrive, for p G (1, oo) and r G [0,1 — 1/p], at 



V-mr,p,E<C\\^\ 



l,p,T 



Vu G VF''P(f ). 



(17) 
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In the case p = 1, we conclude from |Gag57 Thm. l.II] by a similar argument as in |Nec67| Lemma 
2.5.4] that 

II* - ltv\\i-i,i,E < - Itv\\i,i,f Vi) e W''\f), 

hence, 

\\v-lti\l,i,E<CMixf yvGW^'Hf) 
for all r G [0, Z — 1]. The case p = oo is a simple consequence of the fact that the trace operator is 
the classical restriction due to the embedding theorem [Nec67| Thm. 2.3.8]: 



for all r G [0,/]. Thus (17) is proved for all p G [l,oo] with Ip > d and all r G [0,1 — 1/p]. 
Now, if additionally r G No, we have that 

\^-4^\r,p,E ^ II* - ^h*llr.p,£ = inf^ Wild - I^){V + p)\\ E 

pePfc{T) 

< inf J|(/d-/^-)(*+P)ll.,p,E (18) 

< C inf + ^<C7|{;|^ ^. 

peP,_i(T) 

Using the interpolation inequality ( |13[ ) and performing the back-transformation, we get 

< C\\J^%,\det{JE)Mv\i^j^^f 

< c^(ll^Tii,.||j,-MniJTii^r(m)'''H,„T 



where we have used the simple estimate pt Pe together with Q. < 



COROLLARY 2 Let T he an element of an affine partition Th such that the corresponding family 

of partitions {Th}h>o locally quasiuniform. Assume that / G N, p G [l,oo] and Ip > d. Then, for 

jk 



the Lagrange interpolation operator Lf^, there exist constants C > independent of Kt such that 



\v-L^v\r,p,T < Ch'^'''^''^^^'^-'\\v\\i,p^T, 0<r</, (19) 
\v - L^,v\r,p,E < C/^r^'+'''^"'/^-ni^ll/,P,T, 0<r<l-l/p, (20) 

for all V G W^'P{T). 

Proof: For 1 < I < k + 1, the statement coincides with Lemma [3j In the case A: + 1 < /, r G M_|_, 
Lemma [3] implies that 

\V - L^v\r,p,T < Ch^-^^-''\v\k+l,p,T- 



From /i^^^~''|w|fc+i,p,T' < C/i™™^^^"^''^ ''||^||/,p,r the estimate (19) follows. The proof of (20) runs 
analogously. < 

A further possibility of approximating functions in Sobolev spaces is given by the projection w.r.t. 
the L^ inner product. 

DEFINITION 4 The orthogonal L^ -projection : L^{T) Qk{T) is defined, for v G L^{T), 
by 

(v- pI:v,w)o,t = Vu;GQfc(r). (21) 
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In general, the relation (21) is equivalent to a system of linear algebraic equations. It can be 
solved easily provided an L^-orthogonal basis is used. For elements T of an affine partition and 
any multiindex a G Nq, ||a||£oc < k, we have that 

/ dx = I det(JT)| / i'ai'is dx = I det(JT)| TT ^ , . '^a/? = pi^a^p, (22) 

Jt Jf ^Q^i + J- 

where 
and 

denotes the Oj-th one-dimensional Legendre polynomial of degree ai w.r.t. Xj € /, 1 < i < d. In 
this case and with an appropriate indexing, we see that P^u = X^j=i' , ''Ai)o,T'0j- 



Summarized representations about Legendre polynomials can be found in |KS05l App. A] or |QV94 
Ch. 4]. 

From the definition of the L^-projection, the following properties easily follow: 

P'^v = V yvGQkiT), (23) 
\\Ptvh,2,T < \\v\\o,2,T yveL\T). (24) 



LEMMA 5 Let T he an element of an affine partition Th such that the corresponding family of 
partitions {Th}h>o is locally quasiuniform. Then, for 1 < Z < A; + 1, / G N, there exist constants 
C > independent of hx and k such that 

\v-Ptv\r,2,T < C^\v\i,2,T, 0<r<l, (25) 
\v-Ptv\r,2,E < C^pSl7^IH^,2,T, < r < Z - - , < £ « 1 , (26) 



and 



, / + l/2-2r, r>l, 
^("'^)=^^ /-3r/2, 0<r<l, ^^^^ 



for all V e H^''2(r). 



Proof: By |CQ82[ Thm. 2.4], there exists a constant C > independent of k such that, for 

r,/ G M+, r < I, 

\\V - Ptv\\r,2,f ^ C'^~^^''''^ll'^lb,2,f e ^''^(T). 

As in the proof of Lemma [3| this estimate together with (23) and the Deny-Lions lemma ( |Cia91t 
Thm. 14.1]) implies that, for / < A; + 1, 

where C = C{d, I, T) > 0. In particular, the constant does not depend on the polynomial degree k. 
At the faces of T we make use of the following argument. Let E C dT be an arbitrary but fixed face. 
Without loss of generality we may consider it as a subset of M'^^-'^, where the elements x G M.'^^^ 
are characterized by j;^ = and the elements of T satisfy the condition x^ > (otherwise we 
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apply a rotation and a translation of the coordinate system; both operations do not affect the 

differentiability properties of the elements of the function spaces under consideration). 

Now, let w G Vl^*+^/2'2(T) for some s £ (0, / — 1/2]. Then, by Kalyabin's results on Stein's extension 

operator ( [KalSSj p. 185 and Thm. 1]), there exists a total extension operator Ef : W'^~^^^'^''^{T) — )• 
^6+1/2,2 g^pj^ ^^^^ 

\\Efw\\s+l/2,2,R'i < ^II^IL+l/2,2,f- 

The trace theorem ( |AF03| Thm. 7.43 together with Rem. 7.33]) implies that 

l|-^f'^IU2,R<*-i < C\\Efw\\g^i/2,2,R'ij 

SO 

ll-^t^lls,2,_B - ll-^t'^L,2,R'*-l ^ ^ll'"^IL+l/2,2,r- 

Note that for small s G (0,1), we have by a direct trace theorem for Lipschitz domains ( |JK951 
Thm. 3.1]) that 

ll'^IL,2,B - ll'"'IL,2,9f - '^II^IL+l/2,2,f- 

Therefore, Efw\j^ = w\j^ in the sense of W^''^{E) (consider, for I > 2, w £ W^~^^^'^''^{T) for 
s G [1, ^ — 1/2] as an element of, say, w € W^''^{T)) and we finally arrive at the estimate 

\\w\\s,2,E ^ ll^t^L,2,M''-i < C'II^IL+l/2,2,t' « G (0, / - 1/2]. 

In summary, for r E [0, / — 1/2) and < e ^ 1 sufficiently small, we get 

\\w\\r,2,E ^ M\r+e,2,E ^ '^M\r+e+l/2,2,f- 

In particular, for w := v — Pj^v, we have 

11^ - PhV\\r^2,E - ~ ^h^\\r+e+l/2,2,f^ 

and, by |CQ82t Thm. 2.4], 

IP ^h^\\r,2,E - ^'^ ll^ll«,2,T- 



Thus, by (23) and the Deny-Lions lemma ( |Cia9H Thm. 14.1]), 

P ^h^\r,2,E - ^'^ l^lz,2,T- 

The back-transformation to T runs analogously as in the proof of Lemma [3| < 

In contrast to the estimates of the projection error presented above, the next assertion can be 

proved without the use of the Deny-Lions lemma. 

LEMMA 6 Let T he an element of an affine partition Th such that the corresponding family 
of partitions {Th}h>o is locally quasiuniform. Then, for < / and v G W''''^{T) with 1 < s < 
mm{k + 1,1}, there exist constants C > independent of h^ and k such that 

\v-Ptv\o,2,T < C'(x) l''l^'2,T, (28) 
\V - Pl:v\i,2,T < Cp^\vU,2,T, (29) 



V - Pl:v\o^2,dT < C'^^\v\s,2,T- (30) 



-1/2 
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Proof: |(;eon3[ Cor. 3.15, 3.19, (3.1)-(3.3)], |SvdVvDn6[ Lemma 6.1, Rem. 6.2]. < 

REMARK 2 For comparison, we mention the following special case of Lemma^ 

|^-PMo.2,i.<g ^,_3;,(,/,^,) b|/.2,T < C^\v\l,2,T. (31) 

We see that both estimates(30) and (31) are suboptimal of order k^/"^ . The reason lies in the behavior 
of the LP' -projection error measured in the \ -\r^2,T-norm for r > 0. Indeed, any proof of the mentioned 
estimates which is based on a trace estimate will result in right-hand side bounds where the occuring 
norms cause a suboptimal result. 

4 Quadrature and lumping 

In order to evaluate the occuring integrals approximately we will consider interpolatory quadrature 
rules (cf. pOnil Def. 8.1]). 

DEFINITION 5 Let T <Z be a nonempty, compact, connected subset with a Lipschitzian 
boundary (cf lAFOSi Def. A quadrature rule on T with n^^j nodes is characterized 

1. by a set consisting of n\^^ real numbers < cj/, . . . ,a;'^j. > called weights, and 

V ^dof J 

2. by a set Q consisting of nh r points {xi,...,x k } C T, where Xi / Xj if i ^ j, called 

dof 

quadrature nodes. 
The largest natural number k such that 

/ pix)dx = y2u;fp{xi) VpGQfc(r) (32) 

is called the degree of precision or the quadrature order of the quadrature rule. 
From 

p{x)dx = / p{Ft{x))\ det{ JT)\dx 

= j. p{FTixa))<^aix)\det{JT)\dx 

,||a||^oo<fe Q^SNq, ||a||foo <fc 

we immediately get the weights 

= / (^j(x)l det( Jt)| = / Lpi{x)dx, I < i < 
Jf Jt 



dof' 



corresponding to the nodes Xi := Fxixi), where (pi is the Lagrange basis function to Xj. For reasons 
of numerical robustness of the quadrature rule the following condition has to be satisfied: 



w/ > for 1 < i < . (33) 
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The Definition [5] of the nodal quadrature rule can be used to define a discrete inner product. The 
induced discrete norm allows to derive norm equivalence estimates w.r.t. the L^-norm for discrete 
arguments such that the equivalence constants depend only on p, the polynomial degree k and the 
distribution of qudrature nodes. 

To do so, we start with the definition of control volumina and associated lumping operators. The 
discrete L^-norm is given by 



l0,2,T, 



h ■= (.V,v)o,T,h--- 



k k 
ai=0 aa=0 



(34) 



with ^Qgi^d \\a\\icxi<k^a — l^li > 0, Va G Ng. The generalization to the L^-norm for p G [l,oo) 
is straightforward. 

The control volumina are introduced as follows: 



/CKl — 1 



X • • • X 



j=0 



i=0 



^ i=0 i=0 / 



(35) 



where we use the convention ^ 
of Via we have that 

and 



"1 , ,J ._ 



0. As a consequence, for the d-dimensional Jordan measure 



E 



\T\. 



(36) 
(37) 



a&ii. ||a||foo<fc 



The lumping operator ■ C(T) — L°°{T) is defined by 

"dof 



(38) 



i=l 



where 



Xg{x) 



xiG, 
x£G 



is the indicator function of the set G. Integrating the p-th power of the lumping operator, from 
( 36 ) , ( 37 ) together with the affine transformation of the reference element we see that 

'^hivWdx 



\Lhiv)r 



0,p,T 



p dor /* 

|r| / \Lh{v)\Pdx = y^\nj\ / |L, 
Jf ^ Jf 



k 

'dof 



i=0 



"^dof 



dx = |r2j| |t;(3;j)|^ 



(39) 



Furthermore, 



N= 



\0,p,T,h 



Xj)\'^ipj{x) dx 



[ t 

JT 



(40) 
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LEMMA 7 There exist constants Clj^, > independent of Ht such that the following equiv- 
alence estimates are valid: 

ClAv\\o,p,t < \\Lhiv)\\o,p,T < ClM\o,p,t G Qfc(r), p e [2,oo). (41) 

Proof: First we mention that, as a consequence of the substitution rule, the constants do not 
depend on Ht for all elements T which result from an affine transformation of the reference element. 

k k 

Since ||-Z^ft,('y)||opr = J2j=i ll^ll^p = Si=°i obtain, by Nikolski's lemma 

(cf. Lemma 17) and Holder's inequality for sums, 



and 




o,p,r 



\\Lh(.v)\\Q f < max IQj] 

\l<i<nSof / 

< A,nin(M)-V2( 

< A,nin(M)-'/2|T|(P-2)/2P 

where Ajniji(M), A^g^^iM) denote the minimal and maximal eigenvalues, resp., of the corresponding 
mass matrix with entries Mij := Jj, (pi(pj dx. < 

In the approximation of the L^-projection by means of quadrature rules, the choice of the positions 
of the quadrature nodes plays an essential role. 

On the other hand, there is also some freedom in the choice of the node set N in the case of 
Lagrange basis polynomials (^f for < i < A;. In the case of coinciding node sets M and Q we see 
that the definitions of the L^-projection and of the Lagrange basis polynomials result in 

dof 

G = {v-P^v,w)o^T^ = Y.'^^ X'^^^)-^Phv){x^)}w{xi) Vu; G Qfc(r) , (42) 

1=1 

that is v{xi) = (P^w)(xj) for all quadrature nodes. 

k 

Thus the representation {P^v){x) = Ylj=i{'Ph'^)j^j(^) -L^-projection leads to {P^v)i = 

v{xi). Consequently, for the above approximation of the L^-projection, we have the implication 
N=Q^P^ = ll 

The best accuracy can be reached by using the Gauss quadrature rules in the following sense. A 
[k + l)'^-node Gauss quadrature rule yields exact results for polynomials of maximum degree 2k + 1 
(see, e.g., |EG04l Prop. 8.2], |BM97l Sect. 13,14] or | CHQZ07t pp. 448]). The quadrature points 
w.r.t. the domain of integration I = [—1,1] are the zeros of the Legendre polynomial '(/'^^{(x), x G /, 
of degree k + 1. However, the use of Gauss quadrature rules does not lead to optimal estimates of 
the quadrature error w.r.t. the iy''2-norm, as the following result indicates. 

LEMMA 8 For all real numbers r and I, r < I, I > 1, there exists a constant C > independent 
of k such that 

\\v - I^v\\r,2,i < \\v\\i^2,i yv G W''\I), 

where e{r,l) is the function introduced in Lemma^ 
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Proof: See |BM97[ (13.15)]. < 

Alternatively, Gauss-Lobatto quadrature rules can be considered. Here, the quadrature points are 
the zeros of the polynomial 

{x + l){x-l){i^ly{x), xel. (43) 

This quadrature rule is exact for polynomials of maximum degree 2k — 1 (see, e.g., |CHQZ07l pp. 
448]). The boundary points are quadrature nodes. In contrast to the previous lemma, the following 
result is valid. 

LEMMA 9 For all real numbers r and I such that 21 > d+r and < r < 1, there exists a constant 
C > independent of k such that 

\\y - 4L^t,2,f ^ Cf'^'^Mixf ^ W^'\f). (44) 
Proof: See [BM97', Thm. 14.2]. < 

In the application to quadrature, we have the following result (cf. Lemma [7]). 

LEMMA 10 Let Q be the set of Gauss-Lobatto quadrature nodes. Then there exists a constant 
C > independent of hx and k such that the following equivalence estimates are valid: 

\\v\\o,2,T < \\Lh{v)\\o,2,T < C\\v\\o,2,T G Qk{T). (45) 

Proof: First we prove the statement for the reference element. On T, we have 



by (39); the remaining part is a consequence of |CQ82[ (3.9)]. The affine back-transformation to 



the original element shows that the constants are independent of hx- 



REMARK 3 Analogously to Lemma^we can conclude that the above inequalities (45) can be 
extended to the general case p G [2,oo). 

A further interesting property of Gauss-Lobatto quadrature nodes is related with the square sum 
of Lagrange polynomials. 

LEMMA 11 Let ipa{x) be the Lagrange polynomials w.r.t. the Gauss-Lobatto quadrature nodes. 
Then we have the estimate 

^c{xf < 1 Vx G t. (46) 

llall^oo <fc 

Proof: The proof is a consequence of the tensor product representation together with |Fej32[ § 2]: 

k k 

^.(x)2= ^(^^^(xi)2...^^^^(x,)2<L < 

QSNf*, ai=0 «d=0 



|a II nOO <fc 



<1 <1 



5 Projection and interpolation errors w.r.t. Gauss-Lobatto quadrature 
nodes 

The Legendre polynomials used in the representation of the discrete L^-projector possess two im- 
portant properties: On the one hand, they form an orthogonal basis of Qfc(r), and, consequently, 
the corresponding mass matrix is diagonal. On the other hand, there exists a hierarchical decom- 
position of Qfc(r) in the following sense (cf. |EG04l Def. 1.18]). 



13 



DEFINITION 6 (Hierarchical modal basis) A family {Bk}keNoj where Bk denotes a set of 
polynomials, is called hierarchical modal basis if, for all G No, the following properties are 
satisfied: 

1. Bk is a basis o/Q^, 

2. Bk C Bk+i. 

If the L^-projection is discretized by means of Lagrange polynomials for the node sets M = Q, then 
the discrete L^-orthogonality is conserved but, unfortunately, the corresponding Lagrange basis is 
not a hierarchical modal basis. This problem can be resolved by introducing the following notion. 
By N : B ^ J\f we denote that bijective function which assigns a set of polynomials B to the node 
set of the corresponding Lagrange basis. 

DEFINITION 7 (embedded hierarchical nodal basis) A family {BAi^^^^k , where Bj de- 

— ^ — dof 

notes a set of polynomials of maximum degree k, is called embedded hierarchical nodal basis of 
degree K, K gNq, if the following properties are satisfied: 

1. Bk with N{Bk) C N{Bk) is a basis ofQx, 

2. Bk C Bk, N{Bk) = N{Bk), 

3. Bk is a basis ofQk- 
EXAMPLE 1 (d=l) 

(i) The Lagrange polynomials w.r.t. the Gauss-Lobatto quadrature nodes form an embedded hierar- 
chical nodal basis of degree 1 and 2. This follows from the fact that N[Bi) := { — 1, 1} C N{Bk), 
k€n, and N{B2) := {-1, 0, 1} C N{Bk), A; = 2, 4, 6, . . . 

(ii) The Lagrange polynomials w.r.t. the Gauss-Kronrod quadrature nodes {xi}^^-^ =: N{Bk) form 
an embedded hierarchical nodal basis of degree K for k = 2K + 2, K & 'Hq. Given K + 1 Gauss 
quadrature nodes {xjj^Q =: N{Bk), Gauss-Kronrod quadrature rules are defined by adding K -\-2 
nodes such that 

N{Bk) C N{Bk), 

vdx = y^ojivixi) VvGQ3X+4(^) 

(see, e.g., WGGROC^ ). 

As an application we investigate some properties of the so-called fluctuation operator, see the 
definition below. These results are important in the numerical analysis of discontinuos Galerkin 
methods for conservation laws. 

DEFINITION 8 Let {Bj}i^^^ k be an embedded hierarchical nodal basis of degree K, K G Nq. 

— -J — dof 

For Qk(T) = VI^'^{T) ^^(r), where V(^'''{T) := spanBK and Vf^iT) := span{i3fc \ Bk}, the 
projector P^'^ '■ L'^{T) — )• V^'^(T) is defined by 

{v-P^'''v,w)o,T,h = yweV^'\T). (47) 

The operator Pf^{T) := Ld — Pj^'^ in L'^(T), where Id denotes the identity in L'^iT), is called 
fluctuation operator. 
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LEMMA 12 Let k be an embedded hierarchical nodal basis of degree K, K £ No, 

— ^ — dof 

consisting of Lagrange polynomials and suppose Af = Q. Then there is a constant C > independent 
of Ht such that 

\V - Ph'''v\r,2,T < Ch'f'[\v\l,2,T + \PhV\l,2,T], r < I, 

\V - P^'''v\r,2,E < Ch'^'^^-'"[\v\l,2,T + \P^v\l,2,T], 1/2 + r<l 

for all V E W^'P{T). 

Proof: Under the above assumptions we have, for all Lagrange basis polynomials 'f^'^ S ■Sfc, that 

ip^'^{xi) = 6ij, 1 < i,j < Ji-d^f, Xi S N{Bk) = N{Bk)- The definition of the projector P^'^ and the 
property 

w G Vf '^T) =^ w{xi) = 0, Xi G N{Bk) \ N{Bk) (48) 

imply that 



{v - P^' v,w)o^T,h = X^'^i \^v{xi) - {Pf^ ' v){xi)^ w 

i=l 

"dof 

^v{xi) - {P^'''v){Xi)^ w{xi) 



and, thus, 



= T.(Pk''v)^^f'\^^) (Ph'^vh 1 < ^ < nL • (49) 

As a consequence, if denotes the Lagrange interpolation operator w.r.t. the nodes N{Bk), we 
obtain the relation 



1=1 i)J=l i=l 
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It can be used, together with Lemma [3j in the estimation of the projection error as follows: 

I r,K,k I ^1 tK \ , \tK uK,k i 

\V - Pf^' v\r,p,T < \V - Ih v\r,p,T + \Ih V - Pf^ ' v\r,p,T 

^ I tK I I \j-)K,k tK / r>K,k\ i . 



LEMMA 13 Let {'Sj}i< •< fc he an embedded hierarchical nodal basis of degree K, K £ No, 

— ^ — dof 

consisting of Lagrange polynomials and suppose Af = Q. Then: 



(i) -P^''^ is a linear operator, 
(a) WPj^'^vW < C{K,k, II • II) llvllo.oo.T for an arbitrary norm \\ ■ \\ on V^'^{T). 
(Hi) For m G N, f G QkiT), the following estimate is valid: 

{PiiVv),P!,{Vv'^^^-'))o,T,h>0. 
Proof: The assertion (i) easily follows from ( |49[ ): 



P^'''{av + Pw) = ^ [av{xi) + (iw{xi)\ ipf'^ = aPj^'^v + /SP^ 



K.k 

w 



i=l 



Item (ii) results from the estimates 



n 



K 

■dot 



iTjK.k II ^ II / \ K,kii^'sr^i I N||| K,k\ 

\Ph v\\ < 2^ \\v{Xi)(Pi ' II < 2^ I 

i=l i=l 



'^dof 



- X] W'^f'^W ll^llo,oo,T < C{K,k, II • ||)||v||o,oo,T 



i=l 



To verify (iii), we mention the following elementary inequality on T {p := 2 



'm): 



Then 



™5of "fof 



1=1 i=l 

"dof 

^ UJiVv{x^)■VvP~\xi) 

uiVvPl\xi)-VvPl\xi)>Q. 



Interchanging the roles of and ' , we obtain analogously (P/f' (V?;),P;f' (V7;2'"-1))o,t,/i > 

Motivated by Lemma 12 and the error estimate (44) we investigate the local interpolation error of 
the Gauss-Lobatto interpolation operator. 
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LEMMA 14 Let T he an element of an affine partition Th such that the corresponding family of 
partitions {Th}h>o ^•^ locally quasiuniform. Assume that 1 < / < A; + 1, Z G Nq and 21 > d + r. Then, 
for the Lagrange interpolation operator Iq^, there exist constants C > independent of hx and k 
such that 

(h \ 
\v\i,2,T, r<l, (50) 

\v - L^Lv\r,2,E < Ck'i-j) bl/,2,r, < r, 1/2 + r < 1 (51) 

for anO<e<.l and all v £ W^'P{T). 

Proof: Using Lemma [9j the proof runs analogously to the proofs of Lemmata [3] and [5] < 



COROLLARY 3 Let 1 < / < i^' + 1. Under the assumptions of Lemmata\T^ and\14\ there exist 
constants C > independent of hx and k such that 

\V - Pgl^v\r,2,T < C'(^) [\v\l,2,T + \PGLy\l,2,T], < r < 1, 

/, \l-l/2-r 

\V - P^tv\r,2,E < Ck'[-^) [\vk2,T + \P§Lvk2,T], < r, 

/or 1/2 + r < 1, < e < 1 and all v e W^'P{T). 



REMARK 4 As in the proof of Lemma 12, the error estimate will he reduced to an estimate of a 
suitahle seminorm. 

In general, the operators V and iQj^do not commute. An estimate of the commutation error is 
given in the next result. 



LEMMA 15 Under the assumptions of Lemma I4. there exists a constant C > independent of 
hx and k such that 



llcd^v) - VIglv\\o,2,t,gl < C ( ^ ) \v\l,2,T- 



hx 



i-i 



Proof: 



\iI;l{Vv) - VI^lv\\o,2,t,gl ^ C\\lI;l{Vv) - VL^lv\\o,2,t 



< C\\V{v - L^lv)\\o,2,t + C\\Vv - 4l(V^;)||o,2,t 

ffhxV'^ fhxV'^ 
- Hl,2,T + Ci—j \Vv\l-i^2,T 

6 Inverse inequalities 

Inverse inequalities are an important tool in the analysis of finite element methods. Using particular 
properties if the tensor product representation, in this section we derive sharp inverse inequalities 
such that, with exception of the difference of the derivative orders and the spatial dimension. 
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all other information is known explicitly. It is well-known that in the Nikolski inequality the 
dependence of the degree of polynomials cannot be improved, as the example 

P(.)=(^^^)' (52) 

with Tn being the n-th Cebysov polynomial of the first kind demonstrates, see |Tim63l p. 263] 
(although in other papers, e.g. |GNP08t Lemma 2.4] or |CB93l Lemma 1], different statements can 
be found). The representation 

r„(x) = \[{X + - 1)" + {X- \/x2 - 1)" 



shows that 1 — T^{x) has the zeros —1 and 1, thus (52) is really a polynomial. 
Now we present a generalization of Nikolski and Markov inequalities (see |Nik51j . |HST37| Sect. 
Ill]) to the tensor product situation, and we include it into the proof of an inverse inequality given 
in |E(:in4l Lemma 1.138]. 

Throughout this section we restrict ourselves to the case A; > 1. 

LEMMA 16 (local inverse inequality) Given a reference element {T, P, E} such that, for I > 



ne 



0, the embedding P C '°°(T) is satisfied. Let {Th]h<^{o,i\ o, locally quasiuniform family of affi- 
partitions. IfO<m<l, then there exist 

1. for 1 < p, q < oo a constant C = C{1, m,p, q, d, ctq, T, P{T)) > such that 

,, ,, m-Z+d(i-i),, ,, ^ ^ , ^ 

\\v\\l,p,T < Ch^ " '' \\v\\m,q,T yveP{T), (53) 

2. for 1 < q < p < oo and P := Qa,.(T) a constant C = C{1, m,p, d, (Tq) > such that 



The proof will be given after the presentation of the above mentioned Nikolski and Markov inequal- 
ities. 

LEMMA 17 (Nikolski) For < q < p < oo and v G Qfc(T), the following estimate is valid: 

\\<,,,f<{i'i + m~'''''-''''Mo,,f- (55) 

Proof: Due to |DL931 Thm. 4.2.6], in the one-dimensional situation we have the following inequal- 
ity: 

ll^llo,p,/ < {{q + l)A;2)-('/P-'/^) \\v\\o,g,i Vz) G Qfc(/), 
where, as above, / = [—1,1]. Then, for 1 < i < d, we also have that 

11-0(X1, . . . ■,Xi+i, . . . ,Xd)\\o,OD,I 

< {{q + 1)A;2)^^^ \\v{xi, . . . ■,Xi+i, . . .,Xd)\\o,q,i- 

A successive application of this inequality leads to 

Ut'll'^ - = max • • • max Ivfxi, . . . , xw)]"? 

< max-- - max {{q + l)k'^) fj \v{xi, . . . , Xd)\'^dX(i 

< ((g + l)/c2) max - - - max Jjmayi\v{xi, . . . ,Xii)\'^ dxd 



< ((g + 1)A; ) max - - - max fj fj\v{xi, . . . ,Xd)\'^ dxd-idxd 
xiei Xd-2&i 
g 

0,q,f' 



< ■■■ <iiq+l)ky\\v\\'^, 0<q<oo. 
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It remains to make use of an L^'-type interpolation inequality (see, e.g., |EG04l Cor. B.7]): 



— " "0,(?,T ^ ' " "0,g,T 



< {{q + l)k'^) p ll^llogf' 0<q<P<oo. < 



LEMMA 18 (generalized Markov inequality) For v G Q.kil) and p > 1, the following esti- 
mate is valid: 

\\v'\Wi < CMip)k'\\ (56) 

where 

Cm,, = Cm(p) := 2(p - + 1) (i + ^^^f" ■ 

For p = oo, we even have that Cm oo = 1 < lim Cm{p) = 2e. Furthermore, Vp G N : Cm{p) < 
C7m := Ge^+i/^ 

Proof: For 1 < p < oo, the assertion is proved in |HST371 Sect. III]. The estimate Cm{p) < Cm = 
ggi+i/e foj. all p G N can be found in j MMR94l p. 590]. For p = oo, the result follows from |DL931 
Thm. 4.1.4]. < 

REMARK 5 The constants Cm{p) can he improved. For instance, in ]Bar9^ Cor. 2.10] the 
existence of a constant Cm{p) such that lim Cm(p) = 1, P > 2 is demonstrated. 

p— ^-oo 

COROLLARY 4 For 1 < p < oo and v G Qfc(T), the following estimate is valid: 

d I ~ 

Mi^p,f^{ I Y {'^M,pk^)^\\v\\^^p^f (57) 



(with the constant Cm,p as in Lemma 18). 

Proof: Iterating over the order of derivatives, it is sufficient to verify the estimate 

\\d''no,p,f < iCM,pk^)^''^\\v\\o,p,f 

for \a\ = 1. To do so, set a := e^, where denotes the i-th coordinate unit vector. 
For p < oo and 1 < i < d the generalized Markov inequality implies that 

\\v'ixi, . . .,Xi-l, ■,Xi+l,.. ■,Xd)\\lpj 

^ Cli^pk'^^Wvixi, . . . , Xi-i, •, Xi+i, . . . , Xd)\\l^pj, 

hence, on T, 

ll^^-Ollg ^ = Jj . . . Jj\div{x)\P dxidxi . . .dxi-idxi+i . . .dxd 

For p = oo, there is a point y £ T such that 

ll^"^llooof = l'9ii'(y)l = max|aiw(yi, . . . . . . ,yrf)| 

< CM,ook'^\\v{yi, ■ ■ .,yi-i,-,yi+i, ■ ■ ■,yd)\\o,oo,i 

< CM,ook ||^^|lo,oo,t- 



19 



Consequently, 



I " IIP 
\v\ 



n,p,T ^t=u^\a\=in "0,p,T 



<E\=oEM=^{CM,,ki"^'\\v\\: 



0,p,T 



0,p,T 



0,P,T' 



and 



H,oo,T 



max II a w 

0<\a\<l 



lo,oo,t ^ (CA./,c«fc2y«l \\v\\^^^^^ 



< {CM,ook^y\\ 



oo,T' 



Now we are ready to prove the inverse estimates ( 53 ) and ( 54 ) . 
Proof of Lemma [lH 

The first estimate is proved along the lines of the proof of |EG04t Lemma 1.138]. In the particular 
case of a Lagrange reference element we use the same idea of proof but make use of the results 
prepared in this section. 

As usual, first the assertion is proved for m = 0. By (57) and (55), on T we have that 



fd + iy 

i,P,T ^ \ I J 



^ < 



IO,q,T- 



(58) 



In order to get the corresponding estimate for the transformed element T, we make use of ([s]), ([T]) 
and (|8| for < j < / : 



I ^ IP 

\v\ 



d N P 



I " IIP 
\v\ 



d 



'o,g,f ' 



where C{j,p,d,ao) = Cj^d (2\QaoCM,p^ ■ By (4) with / = 0, the back-transformation yields 

1 



\v\j,p,T < "^j C{j,p,d,ao) 



d{---) 



(59) 



Under the assumption Kj} > 1 for < j < /, the corresponding Sobolev norm can be estimated as 

ru ^ I hT\ I hr \ q' [ II IIP 

ll'^llo,<7,T (60) 



< Ei=o{r)[Cis,P,d,ao)[^) 

< r]''')lcij,p,d,ao) 



) I 2(g+l)fc2 ] 



P 



^ n n ' 



|P 

llo,g,r' 



and the assertion is proved for m = 0. Now, let < m < / and a be a multiindex with < |a| < L 
In the case \a\ <l — m, the estimate (|60|) and the relation 



Wv\\Ip,t< E I 

\a\<.l—m 



0,p,T 



|P 

\l-m,p,T 
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imply that 

||'9"u||o,p,T < C{1 -m,p,d,ao) (^^^ ( 2(9+1)^2 ) ll^llo,g,T 

< Cil-m,p,d,ao) (21^)'^^^'^ 

where C{j,p,d,ao) := ('^^j^"') ^^C{j,p,d,ao). This inequahty is also valid for / — m < |a| < /, 
because there exist two further multiindices f3 and 7 such that a = /3 + 7, = I — m and I7I < m. 
As in the first case we get 

Wvh,p,T = \\d^{d^v)\\o,p^T 

< C{1 - m,p,d,ao) \\dMo,,,T 

< Cil-m,p,d,ao)(^) (wMw) ' ' \Mm,q,T. 



2(g+l)fc2 

Since this estimate is valid for < lal < /, the norm llf 7^ can be estimated by using its definition, 
and the statement follows with 

C{l,m,p,d,ao) := Ci-m,d[2vdaoCM,pj I ) \ I ) ' 

The equality sign occurs in the case I = m = and p = q. < 



LEMMA 19 (global inverse inequality) Under the assumptions of Lemma 16, for a locally 
quasiuniform family of affine partitions {Th}he{o,i] ^'^^ for clU v £ Wh, there exist 

1. for 1 < p, q < 00 a constant C = C{1, m,p, q, d, ctq, T, P{T), Cqu) > such that 

( E mLt)' < mi,,t)K (61) 

Ten Ten 

2. for 1 < q < p < CO and P = QkiT) « constant C = C{l,m,p,d,ao) > 

C,^h ^ 

Proof: The proof of the first estimate can be found in |EG041 Cor. 1.141]. With 

C = C{1, m,p, d, (To), 

the second estimate follows from Lemma 16 and from the quasiuniformity of {Th}h>o- 



e/h J^eih 



Taking the p-th root, together with || • Wnv < \\ ■ \\iq the statement follows for p,q ^ 00. The proof 
in the case p = q = 00 immediately follows from the definitions of the norms. 
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